Bound-state third-order optical nonlinearities of embedded germanium nanocrystals 
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. Embedded germanium nanocrystals (NCs) in a silica host matrix are theoretically analyzed to 

' identify their third-order bound-state nonlinearities. A rigorous atomistic pseudopotential approach 

^SJ \ is used for determining the electronic structure and the nonlinear optical susceptibilities. This 

study characterizing the two-photon absorption, nonlinear refractive index, and optical switching 
parameters reveals the full wavelength dependence from static up to the ultraviolet spectrum and 
the size dependence up to a diameter of 3.5 nm. Similar to Si NCs, the intensity-dependent refractive 
index increases with decreasing NC diameter. On the other hand, Ge NCs possess about an order 
of magnitude smaller nonlinear susceptibility compared to Si NCs of the same size. It is observed 
that the two-photon absorption threshold extends beyond the half band-gap value. This enables 
nonlinear refractive index tunability over a much wider wavelength range free from two-photon 
absorption. 

o . 

^ PACS numbers: 42.65.-k, 42.65.Ky, 78.67.Bf 

I. INTRODUCTION 

"j^ . Semiconductor nanocrystals (NCs) benefit from the accumulated knowledge in semiconductor physics and the 

■ maturity of the semiconductor industry as well as the new opportunities provided by the nanoscience, hence they 
I , offer unique optical propertiesFi In particular. Si and Ge NCs attract an increasing attention because of their low- 

''O ' cost and microelectronic-compatible photonic applications ranging from light emitting diodes and lasers to solar cells 
Ch ' and other photonic devicesi^i^ Even though Si and Ge are both Group IV elements, there are a number of notable 
O , differences between them such as the band edge effective mass of the carriers are smaller for Ge, whereas the dielectric 

1 i ' constant is larger which results in bulk exciton radius about five times larger for Ge compared to Si.** As a result, the 

confinement effects will be felt starting from larger sizes. Moreover, Ge is a weakly indirect band-gap semiconductor 
' with the direct to indirect band-gap ratio being 1.2, in contrast to 2.9 in Si. Furthermore, the narrower band-gap of the 
, bulk as well as the NC Ge can be preferred in certain applications to harvest the near infrared part of the spectrum;^ 
Finally, the proximity of direct-gap optical transitions in bulk Ge to the fiber optic communication wavelength of 
\ 1.5 /xm range is particularly important. This has recently stimulated extensive interest; notably tensile-strained Ge 
1 photodetectors on Si platform has been demonstrated^ and a tensile-strained Ge-based laser is proposed.— If the latter 
can also be experimentally demonstrated, this will mark the dawn of the germanium photonics era. 

For all-optical switching and sensor protection applications^ as well as in the absorption of the subband-gap light for 

■ the possible solar cell applicationsy^. the nonlinear refractive index coefficient also known as the optical Kerr index, n2, 
, and two-photon absorption coefficient, /3, are the two crucial third-order optical nonlinearities which play an important 

role. Recent experiments show that Ge NCs have enhanced third-order optical nonlinearitiesj ^°'^^i^^'^'^i^^'^^ However, 
differences in the sample preparation methods, the choice of the matrix, the excitation laser wavelength and the size 
distribution of the NCs contribute to the wide variance within these results as shown in Table I. A number of these 
^ ' investigations have observed two characteristic temporal nonlinear response contributions, distinguished as fast and 
slow, but there is no quantitative agreement among themselves J^ii^ii^ Regarding the origin of the nonlinear response, 
some of these reports have stressed the role of the excited-state contribution produced by the linear absorption^iii^iii 
also the involvement of the trap/defect states was addressed J^ii^ Undoubtedly, more experiments are needed to reach 
to a coherent understanding. On the other hand, to the best of our knowledge, there is no theoretical study identifying 
the wavelength and size dependences of 712 and (3 in Ge NCs. Therefore, a rigorous theoretical work may guide and 
inspire further experimental studies on the foregoing investigations. Moreover, it would help in assessing the potential 
role of Ge NCs, if any, in nonlinear device applications mentioned above. 

In this paper, our aim is to present such a theoretical account concerning 112 and /3 in Ge NCs revealing their 
size scaling and wavelength dependence from static up to ultraviolet region together with a comparison with Si NCs. 
Furthermore, we deal with NCs embedded in a wide band-gap matrix representing silica which is the most common 
choice in the actual structures as can be observed in Table I. Since we do not consider any interface defects, strain 
and thermal effects or the compounding contribution of the excited carriers through linear absorption to the nonlinear 
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processes, our results may serve as a benchmark of the ideal Ge NC bound-state ultrafast third-order nonlinearities. 
In Sec. II we describe the theoretical approach for the electronic structure and the expressions for nonlinear optical 
quantities. The results and discussions are provided in Sec. Ill followed by a brief conclusion. 



TABLE I: The summary of existing experimental studies on the third-order nonlinear optical parameters, n2 and /3 of Ge NCs. 
The sample diameter, D, laser excitation wavelength, Xexc, host matrix and sample preparation information are provided. 
Unspecified data is left as blank. 



Reference D (nm) 


Xexc (nm) 


n2 (cm7GW) P (cm/GW) 


Matrix 


Preparation 


Ref.J^ 


3 


800 


2.7-6.9x10"' 




silica 


ion implantation 


Ref.n. 


6±1.8 


532 


2.6-8.2x10"^ 


190-760 


silica 


cosputtering 


Ref.J^ 


6±1.8 


780 


1.5-8x10"^ 


0.18-0.68 


silica 


cosputtering 


Ref.^S 


5.8, 6.4 


532 




95.4, 143 


alumina E-beam coevaporation 


Ref.^i 




800 


IxlO"'' 


0.04 


silica 


PECVD 


Refl5 


5±2 


820 




1190-1940 


solution 


chemical synthesis 



II. THEORY 



The electronic structure of nanoclusters are accurately obtained routinely by means of density functional theory- 
based ah initio techniquesJ^ However, a several nanometer-diameter NC system including the embedding host matrix 
contains thousands of atoms. This large number currently precludes the use of such ah initio pseudopotential plane 
wave techniques. An alternative route is based on the use of semi-empirical pseudopotential description of the atomic 
environmenl^ii in conjunction with the linear combination of Bloch bands as the expansion basis. ^^■-'^^ In the case of 
embedded Si and Ge NCs, this yields results in well agreement with experimental data for the interband, intraband 
optical absorption'^'^ and the Auger recombination and carrier multiplication.^^ All of these are governed by quantum 
processes taking place over several electronvolt energy range. As a matter of fact, this feature forms an essential 
support for applying the approach to the characterization of the third-order nonlinear susceptibilities up to a photon 
energy of 4 eV. We refer to our previous work for further details on the electronic structure.^- The corresponding 
electronic structure for embedded Ge NCs of different sizes are shown in Fig. [TJ The evolution of the effective gap, 
Eq towards the bulk value (as marked by the gray band) can be observed as the NC diameter increases which is the 
well-known quantum size effect. 

In this work, the electromagnetic interaction Hamiltonian is taken as — er ■ E, in other words, the length gauge 
is used. The third-order optical nonlinearity expressions based on the length gauge have proved to be successful in 
atomic-like systems^ but not in bulk systems because the position operator introduces certain difficulties which can 
actually be overcome.-22. Nevertheless, for bulk systems the velocity gauge has been preferred which on the other hand 
possesses unphysical divergent terms at zero frequency (not present in the length gauge) that poses severe obstacles 
in evaluating the nonlinear optical expressions)^ Hence, we have preferred the length gauge for the evaluation of the 
third-order optical expressions due to resemblance of the band structure of NCs to atomic- like systems (cf. Fig. [T|) . 
The susceptibility expression is obtained through perturbation solution of the density matrix equation of motion^ 
Throughout this work, we distinguish the quantities which refer to unity volume filling factor by an overbar, where 
fv = Vnc/Vsc is the volume filling factor of the NC in the matrix, V^c and Vsc are the volumes of the NC and 
supercell, respectively. The final expression is given by^ 
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where the subscripts {a, 6, c, d} refer to Cartesian indices, uj^ 



+ ujfj + LUa, UJ2 = i-^fj + i^a, ^^1 = are the input 



frequencies, r„m is the matrix element of the position operator between the states n and m, hunm is the difference 
between energies of these states, S is the symmetrization operator,-* indicating that the following expression should be 
averaged over the all possible permutations of the pairs (c, uj^), (6, ujp), and (a, Wq), and finally /„,„ = /„ — /,„ where 



/„ is the occupancy of the state n. The is calculated for m ^ n through r„ 



where toq is the free 



electron mass, and p^m is the momentum matrix element. Hence, after the solution of the electronic structure, the 
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FIG. 1: The energy levels of Ge NCs for different NC sizes. All plots use the same energy reference where the bulk Ge band 
gap is marked by the gray band. 



computational machinery is based on the matrix elements of the standard momentum operator, P, the calculation of 
which trivially reduces to simple summations. 

The above susceptibility expression is evaluated without any approximation taking into account all transitions 
within the 7 eV range. This enables a converged spectrum up to the ultraviolet spectrum. In the case of relatively 
large NCs the number of states falling in this range becomes excessive making the computation quite demanding. 
For instance, for the 3 nm NC the number of valence and conduction states (without the spin degeneracy) become 
3054 and 3314, respectively. As another technical detail, the perfect C^v symmetry of the spherical NCa^S results 
in an energy spectrum with a large number of degenerate states. However, this causes numerical problems in the 
computation of the susceptibility expression given in Eq. ([1]) . This high symmetry problem can be practically removed 
by introducing two widely separated vacancy sites deep inside the matrix. Their sole effect is to introduce a splitting 
of the degenerate states by less than 1 meV. 

When solids are excited with light having a frequency below the band-gap at intensities high enough, third-order 
changes in the refractive index and the absorption are observed due to the virtual excitations of the bound charges. 
Accounting for these effects, the refractive index and the absorption become, respectively, n = uq + n2l, a = ao + (31 
, where uq is the linear refractive index, ao is the linear absorption coefficient, and / is the intensity of the light. n2 
is proportional to Re {x'^^^' } j and is given by^l 



"-2(t^) = r, 2 ^' 



where c is the speed of light. Similarly, /3 is given by=i 

where uj is the angular frequency of the light. Note that Eqs. ([2]) and ([3]) are valid only in the case of negligible 
absorption. The degenerate two-photon absorption cross section ct'^-* (uj) is given by^ 



Jv ll-o^ .J 



E 



5{hWf, - 2hL0), (4) 



where F is the inverse of the lifetime; the corresponding full width energy broadening of 100 meV is used throughout 
this work. The sum over the intermediate states, m, requires all interband and intraband transitions. As we have 
mentioned previously we compute such expressions without any approximation by including all states that contribute 
to the chosen energy window. Finally, a^^\uj) and (3 are related to each other through (3 = 2hujw''^\Lo) . 

Another important factor is the so-called local field effect (LFE) which arise in composite materials of different 
optical properties; the LFEs lead to a correction factor in the third-order nonlinear optical expressions given by;^ 
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FIG. 2: Optical Kerr index at unity filling factor, n2 in Ge NCs as a function of the photon energy for diflerent NC sizes. The 
vertical labels in the ordinates apply to both plots in the same row and the horizontal labels in the lower and upper abscissas 
apply to both plots in the same column. 



where eh and eNc are the dielectric functions of the host matrix and the NC, respectively. 

We fix the local field correction at its static value, since when the correction factor is a function of the wavelength it 
brings about unphysical negative absorption regions at high energies. Further discussions of our model are available 
in our previous worki^ 



III. RESULTS AND DISCUSSION 



We have performed extensive computations on Ge NCs with six different diameters, namely, D — 1.13, 1.47, 1.71, 
2.25, 3, and 3.5 nm. For D — 1.13, and 3.5 nm sizes which correspond to smallest and largest diameters, we have 
calculated the nonlinearities at certain important laser wavelengths. As for the rest, the nonlinearities are computed 
at all frequencies up to 4 eV. For generality, we quote the unity-filling-factor values denoted by an overbar which can 
trivially be converted to any specific realization, however, we should caution that the actual amount of Ge atoms 
forming the NCs is usually a small fraction of the overall excess Ge atoms most of which disperse in the matrix 
without aggregating into a significant NC. In accounting for the LFE the host matrix is assumed to be silica which is 
the most common choice (cf.. Table I). 

Under these conditions, n2 is plotted in Fig. [2] as a function photon energy both in eV (upper abscissas) and in 
units of effective gap, Eq (lower abscissas). As expected, below half Eq (i.e., in the transparency region), there is 
a monotonous behavior and above this value the resonances take over. An intriguing observation is about the sign 
of n2 ■ In general, bulk semiconductors change the sign of n2 above their half Eq values.-^S However, in Fig. [2] we 
observe that in Ge NCs this takes place at even beyond Eq', for the D = 2.25 and 3 nm NCs, the emergence of a new 
resonance is seen to develop which can possibly reproduce this sign change at lower than Eg values for the larger 
NCs. This negative sign of n2 is known to be caused by both the two-photon absorption and the ac Stark effecti^ 
In the size range considered in this work, we believe the former to be more dominant in the negative sign of n2. For 
D = 3 nm NC, and below 1 eV, n2 is of the order of /„ x 10^^ cm^/GW. This value is much larger than the bulk 
valuej^*' When compared to the available Ge NC measurements in Table I, this is in very good agreement with Jie et 
aZiii given the discrepancy in the NC size and the host matrix. On the other hand, other measurements are about 
three orders of magnitude lower for the same quantity. 

The situation is similar in the case of /3 which is plotted against the photon energy in Fig. [3l An important 
observation is that the two-photon absorption onset lies further beyond the half band gap value which is possibly a 
manifestation of the indirect band gap nature of the core medium. We should note that (3 is nonzero (albeit very 
small) down to static values due to band tailing^® which is represented as in our previous work on Si NCs^^ through 
the Lorentzian energy broadening parameter, hT, of 100 meV at full width. For D = 3 nm NC, f3 has a value 
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FIG. 3: Two-photon absorption coefficient at unity filling factor, /3 in Ge NCs as a function of the photon energy for different 
NG sizes. The vertical labels in the ordinates apply to both plots in the same row and the horizontal labels in the lower and 
upper abscissas apply to both plots in the same column. 



10" 



S 10' 

s 

s 
s 

i" 10"' 

0} 



Si NC 



a) 



1.0 



800 nm 



1550 nm 



•M,/ r- 



1.0x10"' 



CO 
CD 



S 1.0x10"'^ 

s" 



— I ' 1 — 

1.5 2.0 



— I — 
2.5 



— I — 
3.0 



IK 

-I- a 
3.5 



1.0x10 - 



NC Diameter (nm) 



b) 



Ge NC 



800 nm 

1550 nm 



1.0 1.5 2.0 2.5 3.0 3.5 
NC Diameter (nm) 



FIG. 4: The size-scaling of the real part of the third-order susceptibility evaluated at two different wavelengths, 1550 nm and 
800 nm, for (a) Si and (b) Ge NCs at unity filling factors. The NCs are embedded in silica matrix. The dashed lines are guide 
to the eye for indicating the overall scaling trend. 



4/^ X 10'^ cm/GW around 2 eV. This value is very high compared to the corresponding bulk value, When typical 
volume filling fraction is taken into account, our values are again in order of magnitude agreement with Gerung et 
al}^ and Jie et alw^ both of which are for somewhat larger NCs. It should be noted that there is an outstanding 
disagreement among the experimental data; for instance, the two most recent experimental data measured at very 
close photon energies differ by five orders of magnitude.— This emerging picture about the large discrepancy on 
the n2 and (3 values calls for further experimental investigations, especially probing the ultrafast response. 

The comparison of the size-scaling trends of the real part of the third-order susceptibility for Si and Ge NCs are 
shown in Fig. [H Two different wavelengths are used, 1550 nm and 800 nm both of which fall below the band gap, 
hence they do not experience any linear absorption. Again unity volume filling factor values are quoted. It can be 
observed that for both Si and Ge NCs, there is a common enhancement trend (dashed lines) as the size is reduced 
especially below 2.5 nm. The oscillations for certain diameters is common to both materials, however, they are more 
pronounced in the Ge NCs. This corroborates with the size scaling of the Auger and carrier multiplication lifetimesi^ 
Another important finding is that third-order susceptibility of Si NCs are more than 20 times larger compared to Ge 
NCs of the same size embedded in the same host matrix. 

For optical switching and modulation applications, one needs large tunability of the refractive index, such as 
through the optical Kerr effect without an appreciable change in the attenuation. Hence, as the figure of merit. 
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FIG. 5: Optical switching parameter, n^lfiX in Ge NCs as a function of the photon energy for different NC sizes. The lower 
and upper abscissas apply to both plots in the same column. 



ri2//3A is proposedi?'^ In our previous work, we have observed that the Si NCs possess much superior figure of merit 
parameters compared to bulk Si.^^ Fig. [5] shows that Ge NCs closely resemble the results of Si NCs.^^ Both Si and Ge 
NCs benefit from the significant blue shift of the onset of the two-photon absorption from the half Eq value, which 
enables the design of such switching or modulation elements over an extended wavelength range. A further indirect 
advantage of this could be the suppression of the optical loss introduced by two-photon absorption generated carriers 
at moderately high pump powers which was a major concern in silicon Raman amplifiers,--- 



IV. CONCLUSIONS 



In summary, we have investigated the wavelength and size dependence of the third-order optical nonlinearities in 
Ge NCs, where our results can serve as a benchmark of the bound-state contribution reflecting the ultrafast response 
of an unstrained perfect sample with no size dispersion. Our computed values for n-i and /3 are in agreement with 
some of the existing experimental data which contain several orders of magnitude disagreement among themselves. 
We observe that below the band gap, there is a common enhancement trend of both the real and imaginary parts of 
the third-order susceptibility as the NC size is reduced. Another important finding is that third-order susceptibility 
of Ge NCs are about an order of magnitude smaller compared to Si NCs of the same size and embedded in the 
same dielectric environment. As in the case of Si NCs, the two-photon absorption threshold extends beyond the half 
band-gap value. This enables nonlinear refractive index tunability over a much wider wavelength range free from 
two-photon absorption. As a final remark, our investigation calls for further experimental work especially to probe 
the ultrafast third-order nonlinear response of Ge NCs. 
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